J!LJ , Ed.: L. Matsson, "Nonlinear Cooperative Phenomena in Biological Systems", 

O^ . Proc. of the Adriatico Research Conference, ICTP, Trieste, Italy, 19-22 August 1997 

(World Scientific, Singapore, 1998) pp. 176-194. 

g ! EFFECTS OF LONG-RANGE DISPERSION IN NONLINEAR 

^ ; DYNAMICS OF DNA MOLECULES 

m ■ 

Yu.B. GAIDIDEF'^ S.F. MINGALEEV 
r^ , " Institute for Theoretical Physics, 252143 Kiev, Ukraine. 









> 

o 



P.L. CHRISTIANSEN, M. JOHANSSON, K.0. RASMUSSEN 
Department of Mathematical Modelling, The Technical University of Denmark, 
(^ ■ DK-2800 Lynghy, Denmark. 

fi , A discrete nonlinear Schrodinger (NLS) model with long-range dispersive interac- 

tions describing the dynamical structure of DNA is proposed. Dispersive inter- 
,„^ actions of two types: the power dependence r~" and the exponential dependence 

*'0' e~^^ on the distance, r, are studied. For s less than some critical value, Scr, and 

similarly for /3 < /3cr there is an interval of bistability where two stable station- 
ary states : narrow, pinned states and broad, mobile states exist at each value 
of the total energy. For cubic nonlinearity the bistability of the solitons occurs 
for dipole-dipole dispersive interaction (s = 3), and for the inverse radius of the 
dispersive interaction /3 < /3cr = 1.67. For increasing degree of nonlinearity, a, the 
critical values Scr and Per increase. The long-distance behavior of the intrinsically 
localized states depends on s. For s > 3 their tails are exponential while for 
_ . 2 < s < 3 they are algebraic. A controlled switching between pinned and mobile 

>^ • states is demonstrated applying a spatially symmetric perturbation in the form 

>P ' of a parametric kick. The mechanism could be important for controlling energy 

(_P ' storage and transport in DNA molecules. 

g\ 

C/3 ' 1 Introduction 

CJ 

C/3 ' Understanding the mechanisms of the functioning of biological macromolecules 

^^, (proteins, DNA, RNA, etc.) remains for now the major challenge in molecular 

biology. One of the most important questions is the understanding of gene 
expression. The expression of a given gene involves two steps: transcription 
and translation. The transcription includes copying the linear genetic infor- 
mation into the messenger ribonucleic acid (niRNA). The information stored 
%\ ' in mRNA is transferred into a sequence of aminoacids using the genetic code. 

mRNA is produced by the enzyme RNA-polymerase (RNAP) which binds to 
the promoter segment of DNA. As a result of the interaction between RNAP 
and promoter of DNA the so-called "bubble" (i.e. a state in which 10-20 base 
pairs are disrupted) is formed. The disruption of 20 base pairs corresponds to 
investing some 100 kcal/mole (0.43 eV)ld. 
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In the framework of a linear model the large- amplitude, motion of the bases 
was supposed to occur due to an interference mechanisma. According to this 
model energetic solvent molecules kick DNA and create elastic waves therein. 
As a result of the interference of two counter propagating elastic waves, the 
base displacements may exceed the elasticity threshold such that DNA under- 
goes a transition ta^ kink form which is more flexible. A similar approach 
was also proposed d'cl. The linear elastic waves in DNA are assumed to be 
strong enough to break a hydrogen bond and thereby facilitate the disruption 
of base pairs. In spite of the attractiveness of this theory which gives at least 
a qualitative interpretation of the experimental dataS there are following fun- 
damental difhculties which to our opinion are inherent to the linear model of 
the DNA dynamics: i) The dispersive properties (the dependence of the group 
velocity on the wave-length) of the vibrational degrees of freedom in DNA will 
cause spreading of the wave packets and therefore smear the interference pat- 
tern. Furthermore, it has been shown eI that the amplitudes of the sugar and 
the base vibrations are rather large even in a crystalline phase of DNA. Since 
the large-amplitude vibrations in the molecules and the molecular complexes 
are usually highly anharmonic their nonlinear properties can not be ignored. 
ii) Molecules and ions which exist in the solution permanently interact with 
DNA. These interactions are usually considered as white noise and their influ- 
ence is modelled by introducing Langevin stochastic forces into the equations 
describing the intramolecular motion. It is well known Q that stochastic forces 
provide relaxation of linear excitations and destroy their coherent properties. 
Equivalently the coherence length (the length of the concerted motions) rapidly 
decreases with increasing temperature, iii) DNA is a complex system which 
has many nearly isoenergetic ground states and may therefore be considered 
as a fluctuating aperiodic system. DNA may have physical characteristics in 
common with quasi-one-dimensional disordered crystals or glasses. However, 
it is known □ that the transmission coefficient for a linear wave propagating 
in disordered chain decreases exponentially with the growth of the distance 
(Anderson localization). In this way it is difficult to explain in the framework 
of linear theory such a phenomenon as an action at distance where concerted 
motion initiated at one end of a biological molecule can be transmitted to its 
other end. 

The above mentioned fundamental problems can be overcome in the frame- 
work of nonlinear models of DNA. Nonlinear interactions can give rise to very 
stable excitations, called solitons, which can travel without being smeared out. 
These excitations are very robust and important in the coherent transfer of 
energy tl. For realistic interatomic potentials the solitary waves are compres- 
sive and supersonic. They propagate without energy loss, and their collisions 



are almost elastic. 

Nonlinear interactions between atomsJp-DNA can give rise to intrinsically 
localized breather-like vibration modes I13U. Localized modes being large- 
amplitude vibrations of a few (2 or 3) particles, can facilitate the disruption of 
base pairs and in this way initiate conformational transitions in DNA. These 
modes can occur as a result of modulational instability of continuum-like non- 
linear modes £3 which is created by energy exchange mechanisms between thje 
nonlinear excitations. The latter favors the growth of the large excitationsOj. 

Nonlinear solitary excitations can maintain their overall shape on long time 
scales even in the presence of the thermal fluctuations. Their robust character 
under the influence of white noise was demonstratedliJ and a simplified model of 
double-stranded DNA was proposed and explored. Quite recently the stability 
of highly localized, breather-like, excitations in discrete uanlinear lattices under 
the influence of thermal fluctuations was investigated ta. It was shown that 
the lifetime of a breather increases with increasing nonlinearity and in this way 
these intrinsically localized modes may provide an excitation energy storage 
even at room temperatures where the environment is seriously fluctuating. 

Several theoretical models have been proposed in the study of the nonlinear 
dynamics and statistical mechanics of DNA (see the very comprehensive reviesi^ 
Ej) . A particularly fruitful model was proposed by Peyrard and Bishop EH 
and Techera, Daemen and Prohofskylla. In the framework of this model the 
DNA molecule is considered to consist of two chains that are transversely 
coupled. Each chain models one of the two polynucleotide strands of the 
DNA molecule. A base is considered to be rigid body connected with its 
opposite partner through the hydrogen-bond potential V(u„), where m„ is 
the stretching of the bond connecting the bases, n, and n — 0,±1,±2,.. is 
labelling the base-pairs. The stretching of the Ti'th base-pair is coupled with 
the stretching of the m'th base-pair through a. dispersive potential K{u„,Um)- 
The process of DNA denaturation was studiedll3'll2l under the assumption that 
the coupling between neighboring base-pairii_is harmonic K{un,Un+i)- An 
entropy-driven denaturation was investigatedcJ taking into account a nonlinear 
potential Khj^^Mn+i) between neighboring base-pairs. The Morse potential 
was chosenllZl'll3'El as the on-site potential V{un) but also the longitudinal wave 
propagation and the denaturation of DNA has been investigated liJ using the 
Lennard-Jones potential to describe the hydrogen bonds. 

In the main part of the previous studies the dispersive interaction K was 
assumed to be short-ranged and a nearest-neighbor approximation was used. 
It is worth noticing, however, that one of two hydrogen bonds which is re- 
sponsible for the interbase coupling: the hydrogen bond in the N — H...O 
group is characterized by a finite dipole moment. Therefore a stretching of 



the base-pair will cause a change of the dipole moment so that the excitation 
transfer in the molecule will be due to transition dipole-dipole interaction with 
a Ijr^ dependence on the distance, r. It is also well knowjx.Lh|at nucleotides in 
DNA are connected by hydrogen-bonded water filamentsE^O. In this case an 
effective long-range excitation transfer may occur due to the nucleotide-water 
coupling. 

In the last few years the importance of the effect of long-range interac- 
tions (LRI) on the properties of nonlinear excitations was demonstrated in 
several different areas of physics. The effective mass of solitons in the Frenkel- 
Kontorova rnndel with a repulsive LRI, their shapes and Emerls barriers were 
investigated a. An implicit form of solitons jsyas obtainedtj in a sine-Gordon 
system with a LRI of the Kac-Baker typecJ'EEl and the dependence of the soli- 
|tan width and energy on the radius of the LRI was analyzed. It was postulated 
EH that the nonlinear term in the sine-Gordon equation has a non-local char- 
acter and novel soliton states, of topological charge zero, were found to exist 
at a large enough radius of the interaction. The effects of long-range inter- 
actions of the Kac-Baker type were studied in static and dynamic nonlinear 
Klein-Gordon modelsEa, and nonlinear Schr6dingerc2l continuum models. The 
effects of a long-range harmonipJnteraction in a chain with short-range anhar- 
monicity were also considered e3. It was demonstrated that the existence of 
two velocity dependent competing length scales leads to two types of solitons 
with characteristically different widths and shapes for two yelocity regions sep- 
arated by a gap. A nonlocal NLS equation was proposed E3 for systems with 
long-range dispersion effects. In contrast to the usual NLS equation station- 
ary solutions only exist for a finite interval of the number of excitations. In 
the upper part of this interval two different kinds of stationary solutions were 
found. The new one containing a cusp soliton was shown to be unstable. It 
was also pointed out that moving solitons radiate with a wavelength propor- 
tional to the velocity. Quite recently Eil we proposed a new nonlocal discrete 
NLS model with a power dependence on the distance of matrix element of 
dispersive interaction. It was found that there is an interval of bistability in 
the NLS models with a long-range dispersive interaction. One of these states 
is a continuum-like soliton and the other is an intrinsically localized mode. 

The goal of this contribution is to investigate the effects of long-range in- 
teractions on the nonlinear dynamics of the two-strand model of DNA. In Sec. 
II we present the analytical theory and the results of numerical simulations 
of stationary states of the discrete NLS model with a long-range dispersive 
interaction. We discuss the bistability phenomenon for the soliton solutions 
and their stability. In the analytical part of this section we use a variational 
approach exploiting an exp-like function as a trial function. Then, in Sec. Ill 



we investigate the long-distance behavior of the nonUncar excitations and show 
that intrinsically localized states of the discrete NLS model with a dispersive 
interaction decaying slower than 1/r^ have algebraic tails. Section IV is de- 
voted to the investigation of switching between bistable states. We show that 
a controlled switching between narrow, pinned states and broad, mobile states 
with only small radiative losses is possible when the stationary states possess 
an internal breathing mode. 

2 System and equations of motion 

We study the two-strand model of DNA which is described by the Lagrangian 

L = T -K~V , (1) 

where 

n ^ ' 

is the kinetic energy (the mass of the base-pair is chosen equal to 1), 

-ft^ = T ^ ^ Jn-raiUra-Unf (3) 

is the dispersive interbase-pair interaction of the stretchings and 

F^^FM (4) 

n 

is the potential energy which describes an intrabase-pair interaction. In Eqs. 
(0)~(q) "- ^^^ "^ ^^^ ^^^^ (base-pair) indices, u„ is the base-pair stretching. 
The value u„ = corresponds to the minimum of the intrabase-pair potential 
V(un). We investigate the model with the following power dependence on the 
distance of the matrix element of the base elastic coupling 

Jn~m = Jl\n-m\\ (5) 

The constant J characterizes the strength of the coupling. The parameter s is 
introduced to cover different physical situations including the nearest-neighbor 
approximation (s — oo), quadrupole-quadrupole (s = 5) and dipole-dipole 
(s — 3) interactions. We shall show that this equation having "tunable" prop- 
erties illuminates both the competition between nonlinearity and dispersion 
and the interplay of long-range interactions and lattice discreteness. To take 



into account the possibility of an indirect coupling between base-pairs (e.g. via 
water filaments) we consider also the case when the matrix element of the base 
elastic coupling has the form 



Jn- 



^ jg-/3|n-m| 



where /3 is the inverse radius of the interaction. 
Assuming that 



d^ViUr,) 



dul 



> 



d'V{u^) 



dul 



for j = 3, 4. 



(6) 



(7) 



i.e. the anharmonicity of the intrabase-pair potential is small, we will use a 
rotating-wave approximation 



■0ne" 



c.c. 



(8) 



where u) 






M„=0 



is the frequency of the harmonic oscillations, ipnit) 



is the complex amplitude which is supposed to vary slowly with time. Inserting 
Eq. @ into Eqs. (|^)-(Q) and averaging with respect to the fast oscillations of 
the frequency uj we conclude that the effective Lagrangian of the system can 
be represented in the form 



-C = 2 51 (^""^n ~ ^nV'n) ~ 'H 



(9) 



Here V'n = H7V'n, T= 2i:- 

n = K. + v 

is the effective Hamiltonian of the system where 



/C = - ^ ^ Jn-,n\i',n - ^n\'^ 



is the effective dispersive energy of the excitation and 



« V 



dtV{^Pne-"^*+C.c)-iO^\lPn? 



(10) 



(11) 



(12) 



is the effective intrabase-pair potential. Usually either a Morse potential P^^H^ 



ppt 



or a Lennard- Jones potential EJ is used to model the hydrogen bonds. With 
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these potentials however it is very compHcated to obtain any analytical results. 
Therefore to gain insight into the problem we will use a simplified nonlinear 
potential in the form 

V--7^El^»l^^^"^ (13) 

where the degree of nonlinearity tr is a parameter which we include to have the 
possibility to tune the nonlinearity as well. 

From the Hamiltonian (nuj we obtain the equation of motion itpn = t: — 

for the wave function i/'„(t) in the form 

#n+ Yl •^n-™(V'm-V'") + IV'nP'^V'n =0 . (14) 

The Hamiltonian Ti and the number of excitations 

^ = Ei^-i' (15) 

n 

are conserved quantities. 

We are interested in stationary solutions of Eq. (14) of the form 

V'n = 4>n exp(iAr) (16) 

with a real shape function 0„ and frequency A. This reduces the governing 
equation for (j>n to 

A0„ = J E l"-'^l"'(0m-0n)+0i'"+'^. (17) 

Thus Eq. (O) is the Euler-Lagrange equation for the problem of extremizing 
Ti under the constraint N = constant. 

To develop a variational approach to the problem we use an ansatz for a 
localized state in the form 



(j)n = VN tanh a exp(— Q:|n|) , (18) 



where a is a trial parameter. The ansatz (h8|) is chosen to satisfy automati- 
cally the normalization condition (n5[) such that the problem of extremizing Ti, 
under the constraint N = constant is reduced to the problem of satisfying the 

■ dn ^ 

equation — — = 0. 
da 
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Figure 1; Number of excitations, N, versus frequency, A, numerically from Eq. (hTl) for 
s = DO (full), 4 (dotted), 3 (short-dashed), 2.5 (long-dashed), 2 (short-long-dashed), 1.9 

(dashed-dotted) . 



Inserting the trial function (18) into the Hamiltonian given by Eqs. (|10D, 
( |ll| ) , (|l^) ^d (^) , and evaluating the discrete sums which enter in these equa- 
tions (see til for details) we get the dispersive part of the Hamiltonian 



/C = 2NJ {C(s) - tanh(a)i^(e~", s - 1) - F(e"", s)} 
and the intrabase-pair potential 



V 



CT+1 



f„ , with /^ = tanh'^+^H coth[((T + l)a] 



In Eq. (19) 



C(.) = ^n- 

n=l 

is the Riemann's zeta function and 

oo 

F(z,s) = ^(z"K) 

n=l 

is the so-called Jonqiere's functionl23. 
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Figure 2: Shows endpoints of the bistability interval for A^ versus dispersion parameter s. 

For s = Scr the endpoints coalesce. Analytical dependence (full), Scr — 2.72. Numerical 

dependence (dashed), Scr — 3.03. 



dn 

da 



According to the variational principle we should satisfy the condition 



which yields 



TV" = 2(cr + l)J(tanh(a) F(e-",s-2) 

+ tanh'(a) F{e-",s-l) )(^ 

\ da 



(23) 



As a direct consequence of Eq. (n7f) , frequency A can be expressed as 



A = -1(/C + 2V) 



(24) 



with /C and V being defined by Eqs. (|l9| ) and (po|). Let us discuss first the 
stationary states of the system for the case a — 1. 

Figure 1 shows the dependence N{A) obtained for a = 1 from direct nu- 
merical solution of Eq. (p7|). A monotonic function is obtained for s > Scr- For 
Scr > s > 2 the dependence becomes non-monotonic (of J\f -type) with a local 
maximum and a local minimum. These extrema coalesce at s = Sc- ~ 3.03. 
For s < 2 the local maximum disappears. The dependence -/V(A) obtained an- 
alytically from Eqs. (E3h and (G4) is in a good qualitative agreement with the 




Figure 3: Hamiltonian of the system, "H, versus frequency, A, numerically from Eq. ( |l7| ) for 

s = 2.5. 

dependence obtained numericallytil. Thus the main features of all discrete NLS 
models with dispersive interaction Jn-m decreasing faster than \n — m\~^'''^ co- 
incide qualitatively with the features obtained in the nearest-neighbor approx- 
imation where only one on-site stationary state exists for any TV. However, for 
2 < s < Scr three stationary states with frequencies Ai(iV) < h.2{N) < A3 (TV) 
there exist for each A'' in the interval \Ni{s)^ 7V„(s)]. In particular, this means 
that in the case of dipole-dipole interaction (s = 3) multiple solutions exist. It 
is noteworthy that similar results are obtained when the dispersive interaction 
is in the form of the Kac-Baker potential (0) . In this case the bistability takes 
place for j3 < 1.67. According to the theorem which was recently provenE3, 
the necessary and sufficient stability criterion for the stationary states is 



dN c! v^ ,2 n 



(25) 



Therefore we can conclude that in the interval [7Vi(s), 7Vu(s)] there are only 

two linearly stable stationary states (Ai(A^) and A3(iV)). The third state is 

1, . dN ^ . . 
unstable smce — — < at A = A2. 
dA 
At the points {A{Ni) and A(iV„)) the stability condition is violated, since 

dN 
(-^-7-)s vanishes. Constructing the locus of the end points we obtain the curve 
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Figure 4: Shapes of the three stationary states for s = 2.5 and N = 3.1. The stable: A : 
(full), A = 0.74 (long-dashed). The unstable: A = 0.57 (short-dashed). 



: 0.21 



presented in Fig. 2. This curve bounds the region of bistabihty. It is analo- 
gous to tha-critical curve in the van der Waals' theory of Uquid-vapour phase 
transitions. Thus in the present case we have a similar phase transition like 
behavior where two phases are the continuum states and the intrinsically local- 
ized states, respectively. The analog of temperature is the parameter s. Figure 
3 shows the multistability phenomenon in terms of the Hamiltonian Ti. of the 
system given by Eqs. (|lO|)-(E2|) for cr = 1 and s — 2.5. For s < s^r there 
is an energy interval where for each Ti. three stationary states with different 
A,£a;ist. The observed bistability is very similar to the recently observed one 
e3C3, where the nearest-neighbor case with an arbitrary degree of nonlinearity 
a was studied. The bistability appears in this case for a above a certain critical 
value. 

Figure 4 shows that the shapes of these solutions differ significantly. The 
low frequency states are wide and continuum like while the high frequency 
solutions represent intrinsically-localized states with a width of several lattice 
spacings. It can be obtained Eil that the inverse widths of these two stable 
states are 



ai 



N 
8J 



l/(s-2) 



TV 
8J 



lnf/(l-21nf) 



as 



M-^ 



(26) 



where (. = exp(l/s) is the characteristic length scale of the dispersive interac- 
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tion which is defined as a distance at which the interaction decreases twice. 
It is seen from these expressions that the existence of two so different sohton 
states for one value of the excitation number, N, is due to the presence of two 
different length scales in the system: the usual scale of the NLS model which 
is related to the competition between nonlinearity and dispersion (expressed 
in terms of the ratio N/J ) and the range of the dispersive interaction £. 

Now we turn to discus stationary states of the discrete NLS model given 
by Eq. (|l7| ) with arbitrary degree of nonlinearity. The main properties of the 
system remain unchanged, but the critical value of the dispersion parameter 
Scr is now a function of a. The results of analytical consideration confirmed 
by simulation show that Scr increases with increasing a. In particular, for 
cr > 1.4 (the value at which discrete sjonmetric ground state can be unstable in 
the nearest-neighbor approximationEj) the bistability in the nonlinear energy 
spectrum occurs even for s < 6. 

3 Tails of intrinsically localized states 

Investigating the asymptotic behavior of the excitations, it is convenient to 
rewrite Eq. ( p7| ) (we consider here the case cr = 1) in the form 

</-« = 5]G„-,„(A)(/.3 , (27) 

where 

^ / . ^ 1 f ,, cos(kn) 

is the Green's function with the spectrum function 

oo 

C{k) = 2 ^ J„(l - cos(fcn)) . (29) 



Deriving the asymptotic expressions for the Green's function (^8|)Eil we obtain 
that the tails of the intrinsically localized states are given by the expressions 




'<A + l«W„„f_ji^H . .>3, (30, 



2 < s < 3 , (31) 
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Figure 5: Amplitude in the tail of the stationary state for , 

(circles) and Eq. ^ (full). 



2.5 and n = 450. Numerical 



for \n\ — > cx). Thus we can conclude here that only in the case of the short- 
range dispersion (s > 3) the tails of intrinsically localized states have a usual 
exponential form. In the systems with long-range dispersive interactions these 
states have algebraic tails. Figure 5 shows the long-distance behavior of the 
intrinsically localized states for s — 2.5 and different values of the frequency A. 
It is seen that the form of the tails predicted by Eq. (^) is in a good agreement 
with the results of numerical simulations. 

The long-distance behavior of the intrinsically localized states may play an 
essential role in the thermodynamic properties of the DNA molecule because 
in systems where the interaction decays algebraically it can be responsible for 
the appearance of new thermodynamically stable states (for example, quite 
recently EJ the existence of the Neel order in the ground state of Heisenberg 
antiferromagnetic chains with algebraic long-range interactions was proven). 

4 Switching between bistable states 

Having established the existence of bistable stationary states in the nonlocal 
discrete NLS system, a natural question that arises concerns the role of these 
states in the full dynamics of the model. In particular, it is of interest to 
investigate the possibility of switching between the stable states under the in- 
fluence of external perturbations, and what type of perturbations could be used 
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to control the switching. Switching of this type is important in the description 
of nonhnear transport and storage of energy in biomolecules hke the DNA, 
since a mobile continuum-like excitation can provide action at distance while 
the switching to a discrete, pinned statejpan facilitate the structural changes of 
the DNAE^IIj. As it was shown recentlyE3 the switching will occur if the system 
is perturbed in a way so that an internal, spatially localized and symmetric 
mode ('breathing mode') of the stationary state is excited above a threshold 
value. 

To investigate the time evolution of an initially small perturbation e„(0) 
of the stationary state ( Jl^ ) we write 

V-nW = (</)„ + £„ (r))e'A^ (32) 

Decomposing e„(T) into real and imaginary parts, €„ and e„ , we obtain from 
Eq. ( |l4| ) with cr = 1 in the linear approximation 

=-^(f )- ( _v ;■ 11 1) 1 (33) 

where, for a system with M sites, i/+ and H~ are M x M matrices defined 

by 

//± - (A - (2 T l)(t>l + 2) S,,, - J,_„ (34) 

with Jo = 0. By definition the stationary solution is linearly stable if the 
perturbation e„(T) as calculated from Eq. ( p3[ ) remains bounded for all times. 
Linear stability is then equivalent to the matrix A4 having no eigenvalues with 
a positive real part. Changing some parameter (e.g. A), a stable state might 
become unstable. The 'direction' in which an initial perturbation will grow 
is then determined by the eigenvector corresponding to the eigenvalue of M 
with a positive real part. We will in sequel mainly discuss the case when the 
matrix element of base elastic coupling Jn-m decreases exponentially with the 
distance |n — to| (see Eq. (g)) with the inverse radius of the interaction [3 — 1. 
For such value of /3 the multistability occurs in the interval 3.23 < N < 3.78. 
It is worth noticing, however, that the scenario of switching described below 
remains qualitatively unchanged for all values of /3 < 1.67, and also for the 
algebraically decaying dispersive coupling with 2 < s < 3.03. 

The study lJ of the eigenvalue problem for the matrix Ai showed the 
existence of a spatially symmetric internal breathing mode for both the nar- 
row and broad components of the bistable state. Furthermore, the low fre- 
quency (broad) component also possesses a spatially antisymmetric transla- 
tional (" pinning" ) mode Ea. Since the appearance of a translational mode 
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Figure 6: Switching from continuum- like to discrete state for (3 = 1. The initial state (f>„ 

has the frequency A ~ 0.31 and N = 3.6. The time evolution of |V'no(''")l^ when a phase 

torsion is applied to the center site with 8 = 0.261 (lower curve) and 9 = 0.262 (upper curve), 

respectively; inset shows time evolution of |'/'n(''")P for 9 = 0.262. 



implies that the stationary state gains mobihty l3, the continuum- hke state 
wiU have a high mobihty. 

An iUustration of how the presence of an internal breathing mode can af- 
fect the dynamics of a slightly perturbed stable stationary state is given in Fig. 
6. To excite the breathing mode we apply a spatially symmetric, localized per- 
turbation, which we choose to conserve the number of excitations in order not 
to change the effective nonlinearity of the system. The simplest choice, which 
we have used in the simulations shown here, is to kick the central site no of the 
system at t = by adding a parametric force term of the form 06n,no^(j')4'n{T) 
to the left-hand-side of Eq. (jlj) . A possible physical motivation of the appear- 
ance of such kind of parametric kick may be the following. It is well known 
that biomolecules in cells interact with solvent molecules and ions (ligands). 
At the position where the ligand links to a biomolecule the potential energy of 
the biomolecule changes. This manifests itself in a local change of the vibration 
frequency of biomolecular units. In our case it means that the frequency tu of 
the internal base-pairs oscillations locally changes and the potential energy (0) 
should be replaced by 



1 



y-^ yK) + -<5„„„0(rX 



(35) 
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where uq is the site where the ligand interacts with the biomolcculc, is a 
constant which characterizes the energy interaction between the hgand and 
biomolecule. It may be a function of time because the hgand may attach or 
detach at different time moments. If we use the rotating wave approximation 
(^ the NLS model given by Eq. (g) with the additional term in Eq. (p^ ) being 



Assuming that 



V„t=^<5„„„^(T)|V'n(r)p. (36) 



W=E^^'^(^-^^-) (37) 



we obtain the model where the attachment and detachment of ligands are 
considered as kicks which occur at the time moments tj. 

As can be easily shown, this perturbation affects only the site np at r = 0, 
and results in a 'twist' of the stationary state at this state with an angle 6, 
i.e. V'no(O) = 4>no 6**. The immediate consequence of this kick is, as can been 
deduced from the form of Eq. (f#, that -^ (IV'noP) '^i^ be positive (negative) 
when 9 > { 9 < 0). Thus, to obtain switching from the continuum- like 
state to the discrete state we choose 9 > 0, while we choose 9 < when 
investigating switching in the opposite direction. We find that in a large part 
of the multistability regime there is a well-defined threshold value 9th, such 
that when the initial phase torsion is smaller than 9th, periodic, slowly decay- 
ing 'breather' oscillations around the initial state will occur, while for strong 
enough kicks (phase torsions larger than 9th) the state switches into the other 
stable stationary state. 

It is worth remarking that the particular choice of perturbation is not 
important for the qualitative features of the switching, as long as there is a 
substantial overlap between the perturbation and the internal breathing mode. 
We believe also that the mechanism for switching described here can be applied 
for any multistable system whereJJie instability is connected with a breathing 
mode. For example, we observedEj a similar switching behavior in the nearest 
neighbor iliscrete NLS equation with a higher degree of nonlinearity a, which 
is known LJ to exhibit multistability. 

5 Conclusion 

We have proposed a new nonlocal discrete nonlinear Schrodinger model for 
the dynamical structure of DNA with long-range (r"" and e~^^) dispersive 
interaction. We have shown that there is a multistability in the spectrum of 
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stationary states of the model with a long-range dispersive interaction s < Scr 
(/3 < (3cr)- There is an energy interval where two stable stationary states exist 
at each value of the Hamiltonian H. One of these states is a continuum-like 
soliton and the other one is an intrinsically localized mode. The existence of 
the bistability phenomenon in the NLS models with a nonlocal dispersion is a 
result of the competition of two length scales which exist in the system: the 
scale related to the competition between nonlinearity and dispersion, and the 
scale related to the dispersion interaction. 

We found that the critical value of the dispersion parameter Scr for the 
on-site stationary state in the case of cubic nonlinearity exceeds 3. This means 
that the bistable behavior may occur in the case of DNA where the stretching 
motion of base-pairs is accompanied by a change of their dipole moments. 

We have shown that the long-distance behavior of intrinsically localized 
states in discrete NLS models with a nonlocal dispersion depends drastically 
on the value of the dispersive parameter s. Only for short-range dispersions 
the excitation wave functions decay exponentially. In the systems where the 
matrix element of base elastic coupling depends on the distance slower than 
1/r^ the nonlinear excitations have algebraic tails. The long-distance behavior 
may be important for the thermodynamics of DNA since it provides long-range 
order in one-dimensional systems. 

We have shown that a controlled switching between narrow, pinned states 
and broad, mobile states is possible. Applying a perturbation in the form of 
parametric kick, we showed that switching occurs beyond some well-defined 
threshold value of the kick strength. 

The particular choice of perturbation is not important for the qualitative 
features of the switching, as long as there is a substantial overlap between 
the perturbation and the internal breathing mode. Thus, we believe that the 
mechanism for switching described here can be applied for any multistable sys- 
tem where the instability is connected with a breathing mode. The switching 
phenomenon could be important for controlling energy storage and transport 
in DNA molecules. 
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